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Abstract 

The energy flow created in pp collisions at y/s = 7 TeV is studied within the pseu- 
dorapidity range 1.9 < rj < 4.9 with data collected by the LHCb experiment. The 
measurements are performed for inclusive minimurmbias interactions, hard scatter- 
T — I ing processes and events with an enhanced or suppressed diffractive contribution. 

The results are compared to predictions given by PYTHiA-based and cosmic-ray 
event generators, which provide different models of soft hadronic interactions. 
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1 Introduction 



In Quantum Chromodynamics (QCD), the final state of an inelastic hadron-hadron col- 
lision can be described by contributions from hard and soft scattering occurring between 
the constituents of the hadrons, initial- and final-state (gluon) radiation and the fragmen- 
tation of the initially coloured partonic final state into colour-neutral hadrons. The soft 
component of a collision is called the underlying event. Its precise theoretical description 
remains a challenge, while the dynamics of hard scattering processes is well described by 
perturbative QCD. One source of the underlying event activity is multi-parton interac- 
tions (MPI). These arise mainly in the region of a very low parton momentum fraction, 
where parton densities are high so that the probability of more than a single parton-parton 
interaction per hadron-hadron collision is large. MPI effects become increasingly impor- 
tant at LHC collision energies, where inelastic interactions between very soft partons are 
sufficiently energetic to contribute to final state particle production [I]. 

MPI phenomena can be probed by measuring in the centre-of-mass system the energy 
flow created in inelastic hadron-hadron interactions at large values of the pseudorapidity 
T) = — ln[tan(#/2)], with 9 being the polar angle of particles with respect to the beam 
axis. In this region of phase space, the energy flow is expected to be directly sensitive 
to the amount of parton radiation and MPI [2]. For a particular pseudorapidity interval 
with width Ar], the total energy flow d-Etotai/ drj, which is normalised to the number of 
inelastic pp interactions N int , is defined as 



1 dE tota \ 
N- mt drj 




(1) 



where N paxtjV is the total number of stable particles and E iiV is the energy of the individual 
particles. 

In this study, the energy flow is measured in pp collisions at ^/s = 7TeV within 
the pseudorapidity range 1.9 < rj < 4.9. This extends the previous measurements that 
have been made in pp [3] and ep collisions (2) to larger pseudorapidity values and higher 
centre-of-mass energies, and complements the study performed by the CMS collabora- 
tion |5|. Experimental results are compared to predictions given by PYTHlA-based [6}[7] 
and cosmic-ray event generators |8j[9], which model the underlying event activity in dif- 
ferent ways. In order to probe various aspects of multi-particle production in high-energy 
hadron-hadron collisions, the measurements are performed for the following four classes of 
events: inclusive minimum-bias, hard scattering, diffractive, and non-diffractive enriched 
interactions. 



2 The LHCb detector 

The LHCb detector [To] is a single-arm forward spectrometer with an angular coverage 
from 10 mrad to 300 (250) mrad in the bending (non-bending) plane, designed for the 



1 



study of b- and c-hadrons. The detector includes a high precision tracking system con- 
sisting of a silicon-strip vertex detector (VELO) surrounding the pp interaction region, 
a large-area silicon-strip detector located upstream of a dipole magnet with a bending 
power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes 
placed downstream. The VELO has a larger angular acceptance than the rest of the 
spectrometer, including partial coverage of the backward region. It allows reconstruction 
of charged particle tracks in the pseudorapidity ranges 1.5 < 77 < 5.0 and —4 < r\ < —1.5. 
The combined tracking system has a momentum resolution Ap/p that varies from 0.4% at 
5 GeV/ c to 0.6% at 100 GeV/ c, and an impact parameter resolution of 20 \xm for tracks with 
high transverse momentum, p^. Charged hadrons are identified using two ring-imaging 
Cherenkov detectors. Photon, electron and hadron candidates are distinguished by a 
calorimeter system consisting of scintillating-pad and preshower detectors, an electromag- 
netic calorimeter (ECAL) and a hadronic calorimeter (HCAL). The calorimeters have an 
energy resolution of a(E)/E = 10%/ VE © 1% and a(E)/E = 69%/ VE © 9% (with E in 
GeV), respectively. Muons are identified by a system composed of alternating layers of 
iron and multiwire proportional chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies a full event reconstruction. 
For the minimum-bias data used in this analysis, the hardware trigger was accepting all 
beam-beam crossings, while the presence of at least one reconstructed track was required 
in the software stage to record an event. 



3 Data analysis 

3.1 Data and Monte Carlo samples 

The analysis is performed using a sample of minimum-bias data collected in pp collisions 
at y/s = 7TeV during the initial running period of the LHC with low interaction rate. 
The fraction of bunch crossings with two or more collisions ( "pile-up events" ) is estimated 
to be approximately 5%. The total number of events available in the sample is 5.8 x 10 6 , 
corresponding to an integrated luminosity of about O.lnb -1 . 

Fully simulated minimum-bias pp events at = 7 TeV were generated using the LHCb 
tune 11 of the Pythia 6.4 event generator (6). Here, decays of hadronic particles are 



described by EvtGen 12 in which final state radiation is generated using Photos 13 



The interaction of the generated particles with the detector and its response are imple- 



mented using the Geant4 toolkit 14,15 as described in Ref. 16 . Additional Monte 
Carlo (MC) samples with fully simulated minimum-bias pp interactions at -y/s = 7 TeV 
were generated using the Perugia and Perugia NOCR [17] tunes of Pythia 6.4. These 
models along with the LHCb tune use different values for the MPI energy scaling parame- 
ter and MPI p^ cut-off, which entails a sizeable deviation in the amount of MPI predicted 
by these tunes. 

The LHCb tune utilises the CTEQ6L parton density functions (PDFs) [18], while 
both Perugia tunes use the CTEQ5L PDFs |19|. Colour reconnection effects are not 
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included in the Perugia NOCR tune. In the MC samples generated with the Perugia 
and Perugia NOCR tunes, diffractive pp interactions are suppressed at generator level, 
whereas the sample generated with the LHCb tune contains the contributions from both 
single and double diffractive processes. A sample of fully simulated diffractive events 
generated with Pythia 8.130 [7], which utilises the CTEQ5L PDFs, is used in addition. 
This event generator gives a more accurate description of diffractive pp interactions than 
Pythia 6, especially at high-px, as it includes the contribution from hard diffractive 
processes, which is absent in Pythia 6 |20| . 

In addition to the models above, experimental results are compared to generator level 
predictions given by the Pythia 8.135 model with default parameters. Furthermore, the 
measurements are compared with predictions given by the cosmic-ray interaction models 



Epos 1.99 21 , QgsjetOI, QgsjetII-03 22 , and Sibyll 2.1 [23], which are widely 



used in extensive air shower simulations and are not tuned to LHC data. These gen- 
erate inelastic pp interactions taking into account the contributions from both soft and 
hard scattering processes. Soft contributions are described with Gribov's Reggeon field 
theory [24] via exchanges between virtual quasi-particle states (Pomerons), while hard pro- 
cesses are described by perturbative QCD via exchanges of hard or semi-hard Pomerons. 
The predictions given by these models diverge mainly because of different treatments of 



non-linear interaction effects related to parton saturation [25] and shadowing 26 . The 
QgsjetOI model describes hadronic multiple scattering processes as multiple exchanges 
of Pomerons without specific treatment of saturation effects. A distinct feature of the 
QgsjetII model is the treatment of non-linear parton effects via Pomeron interactions 
taking into account all order re-summation of the corresponding Reggeon field theory di- 



agrams. Based on the dual parton model 27 , Sibyll utilises the Lund string model 28 



for hadronisation and describes soft and hard processes using the Pomeron formalism and 



the minijet model 29 , correspondingly. The treatment of non-linear effects in this model 
is based on a simple geometrical approach of parton saturation. The Epos model takes 
into account energy-momentum correlations between multiple re-scatterings and describes 
non-linear effects using an effective treatment of lowest order Pomeron-Pomeron interac- 
tion graphs. It also accounts for the final state interaction of the produced particles. 



3.2 Analysis strategy 

The energy flow, as defined in Eq.[TJ is the energy-weighted pseudorapidity distribution 
of particles, normalised to the number of inelastic interactions and the 77-bin size. The 
measurements are performed in ten equidistant pseudorapidity bins of width At] = 0.3 
over the range 1.9 < rj < 4.9. The primary measurement is the energy flow carried 
by charged particles (charged energy flow). It is performed with reconstructed tracks 
which contain hits in the VELO and downstream tracking stations and have momentum 
in the range 2 < p < lOOOGeV/c. Particle identification is not required in this analysis, 
as the energy is taken from the momentum, neglecting particle masses. In order to be 
able to compare the results of the measurements with generator level predictions, the 
reconstructed charged energy flow is corrected for detector effects. The total energy flow 
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is determined by using a data-constrained MC estimate of the neutral component based 
on information from the ECAL, while the HCAL is not used. Details of the procedure 
are discussed below. 



3.3 Event classes 

The event classes studied in this analysis are defined as follows. Inclusive minimum-bias 
events are selected by requesting the presence of at least one track originating from the 
luminous region in order to suppress pollution from beam-gas interactions and beam halo 
related background. Events with two or more reconstructed primary vertices are rejected 
to suppress pile-up contamination. To minimise biases on the track multiplicity of the 
event, the information on the primary vertex is not used. The selected inclusive minimum- 
bias interactions are further classified as hard scattering, diffractive and non-diffractive 
enriched events using the following criteria: 

• Hard scattering events: at least one track with p^ > 3 GeV/c and 1.9 < rj < 4.9. 

• Diffractive enriched events: no tracks reconstructed with —3.5 < i] < —1.5. 

• Non-diffractive enriched events: at least one track reconstructed with —3.5 < rj < 
-1.5. 

The selection requirements applied for the last two event classes are motivated by the fact 



that a sizeable rapidity gap is an experimental signature of diffractive processes 30 . The 
level of enrichment of the diffractive and non-diffractive samples was studied in simulation, 
by retrieving the Pythia process type of the pp interaction for every selected diffractive 
and non-diffractive candidate. In the case of the LHCb tune of Pythia 6.4, the purities 
of the selected diffractive and non-diffractive enriched samples are found to be about 70% 
and 90%, respectively. Although the actual percentages are only meaningful within the 
specific model, the study shows that the applied selection criteria indeed lead to sizeable 
enhancement of the respective event classes. 

To minimise the experimental corrections, the definition of the event classes at gener- 
ator level is similar to that at detector level. Inclusive minimum-bias events at generator 
level are selected by requiring the presence of at least one outgoing final-state charged 
particle (lifetime r > 10 -8 s) in the pseudorapidity range 1.9 < rj < 4.9, but without 
imposing any condition on its energy. The sample of hard scattering events is selected by 
requesting at least one final-state charged particle with p T > 3 GeV/c and 1.9 < rj < 4.9. 
The absence or presence of at least one final-state charged particle in —3.5 < 77 < — 1.5 is 
used as criterion to select diffractive and non-diffractive enriched events among inclusive 
minimum-bias interactions, respectively. For the selected events, the energy flow at gener- 
ator level is determined using the outgoing final-state charged and neutral particle^] which 
are either prompt, originating directly from the fragmentation, or the decay products of 
unstable particles. Since neutrinos are not reconstructed by the LHCb spectrometer the 



"These include ir^ 1 , K ± , e ± , fj^ , p, p, 7, n, n and 
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energy carried by these particles is not taken into account. Only MC events simulated 
with exactly one inelastic pp interaction are considered in this study. 



3.4 Corrections 

The reconstructed charged energy flow measured with data is corrected for detector effects 
using bin-by-bin correction factors, which are estimated as the ratio of the charged energy 
flow at generator and detector level in simulation for each rj region and event class under 
consideration. The overall correction factor for each bin is taken as the average of the 
correction factors obtained with different MC models used in this analysis. For inclusive, 
hard scattering and non-diffractive enriched events, the average and standard deviation of 
the correction factors, which is included in the model-dependent systematic uncertainty, 
are determined from the LHCb, Perugia and Perugia NOCR tunes of Pythia 6.4. In 
the case of the diffractive enriched event class, only the LHCb tune and the Pythia 8 
diffractive simulation are used. Except for the lowest rj bin, which suffers from reduced 
acceptance for low-px particles and thus exhibits large corrections and a sizeable model 
dependence, the correction factors are found to be stable among the models with a slight 
rise towards the edges of the detector acceptance. The majority of the factors are well 
below two, indicating that most of the energy is measured by the detector. In the case 
of diffractive enriched events, the correction factors obtained with the LHCb tune are 
slightly smaller than unity for some of the bins, i.e. the energy flow at detector level is 
found to be larger than at generator level. This is due to detection inefficiency for charged 
particles over the pseudorapidity range —3.5 < r] < — 1.5. As a result, some of the events 
containing backward going charged particles migrate into the diffractive sample, which 
leads to enhanced energy flow at detector level. 

For the measurement of the total energy flow, the neutral component F neut ^ v is esti- 
mated in the following way. To first order F neut r? is assumed to be proportional to the 
corrected charged energy flow -F char r? with a factor R gen ,n, which is the average ratio of 
the neutral energy flow to the charged energy flow obtained at generator level for each r\ 
bin and event class. This ratio is found to be rather stable over the entire pseudorapidity 
range of the measurements with only small variations between the Pythia tunes. This 
reflects the usage of the same hadronisation mechanism governed in the Pythia genera- 
tor by the Lund string model |6j[7]. The latter successfully describes the hadronisation of 
quarks and gluons emerging from high energy interactions and was rigorously tested for 



high-px processes 31 33 . The -Rgen.r? ratio is found to be around 0.6 for all event types 
except the hard scattering interactions. For the latter, it is about 15% smaller for all r\ 
bins. This feature is found to be a consequence of the requirement of a high-pp charged 
particle in the definition of this event class. 

Under the assumption outlined above, the total energy flow for a particular event class 
and pseudorapidity bin F tota i jr? can be written as 

-^total,r/ Fchai,ri -Fneut,^ -^char,^ X (1 -)- -Rgen,r?) • (2) 

In order to constrain this initially purely model-based estimate of the neutral energy flow 
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to data, the total energy flow is further multiplied by an additional correction factor k v . 
It accounts for differences between simulation and data being defined for every r\ bin as 



7 1 ~l~ -^data,ri / \ 

k *> = TTr — • (3) 

Here, -Rdata,r) and -R mcfl are ratios of the uncorrected neutral to charged energy flow mea- 
sured in data and simulation, respectively. The R mc>v ratio is obtained with different 
Pythia tunes and its average is taken for the estimation of the k v factors. The neutral 
component of -Rdata,r) and R mc ^ is measured using reconstructed photon candidates which 
are selected from neutral clusters in the ECAL with an energy greater than 2 GeV and 
Pt > 0.2GeV/c. Since the polar angular coverage of the ECAL begins at about 30mrad, 
there are no measurements of the neutral energy for the last two rj bins (rj > 4.3). The k v 
factors for this pseudorapidity region are estimated using a linear extrapolation of the k v 
factors obtained for the pseudorapidity interval 3.1 < rj < 4.3. The bins with rj < 3.1 are 
not considered for the extrapolation, as these are affected by the detection inefficiency for 
low-px charged particles. The majority of the k v factors are found to be close to unity, 
reflecting the fact that the ratio of the neutral to charged energy flow is well simulated at 
detector level. 



4 Systematic uncertainties 

The total uncertainties on the results are dominated by systematic effects, as the statistical 
uncertainties are found to be negligible for all rj bins and event classes. The various 
contributions to the systematic uncertainties are summarised in Table [TJ 

For all event types except hard scattering interactions, the largest uncertainty on the 
charged energy flow arises from the model dependence of the bin-by-bin correction factors. 
Here, the largest impact is at low rj, reaching 9% for inclusive and non-diffractive enriched 
events, 4% for hard scattering interactions and up to 43% for diffractive enriched events. 
At large rj this effect generally drops to about 1-2% for all event classes except diffractive 
enriched interactions for which it stays above 15%. 

Systematic uncertainties related to the track selection requirements are estimated by 
comparing the fraction of the energy flow from tracks which are rejected by the selection 
cuts in data and simulation. For the majority of the bins the resulting systematic un- 
certainty is found to be less than 4%. Only for hard scattering events this uncertainty 
approaches 9% at low rj. 

To account for differences between the true tracking efficiency and that estimated 
using simulation, a global 3% systematic uncertainty is assigned across the entire rj range 



following the analysis presented in Ref. 34 . This applies for all event classes under 
consideration. 

The other tracking related factors having an influence on the charged energy flow 
measurements are contaminations from multiply reconstructed tracks and tracks created 
from random combinations of hits (spurious tracks). The impact of the former is estimated 
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Table 1: Relative systematic uncertainties (in percent) affecting the energy flow measurements 
for all event classes. The total uncertainties are obtained by adding the individual sources in 
quadrature. The ranges indicate the variation of the uncertainty as a function of rj. 



Source of 


Inclusive 


Hard 


Diffractive 


Non- diffractive 


uncertainty 


minbias 


scattering 


enriched 


enriched 


Model uncertainty on 


0.6-9.2 


0.7-4.1 


16-43 


0.7-8.6 


correction factors 










Selection cuts 


1.0-4.9 


2.7-8.8 


0.9-2.8 


1.1-5.0 


Tracking efficiency 


3 


3 


3 


3 


Multiple tracks 


1 


1 


1 


1 


Spurious tracks 


0.3 - 1.2 


0.4-1.7 


0.2-0.7 


0.3-1.2 


Magnet polarity 






2.6-7.7 




Residual pile-up 


1.7 


1.7 


1.7 


1.7 


Total on F chaXjV 


3.9 - 11 


4.9 - 10 


16-43 


4.0 - 11 


Variation of R gcniV 


0.8-6.1 


0.7-2.9 


1.5-23 


0.9-5.5 


and k v factors 










Photon efficiency 


1.4-1.6 


1.2 - 1.3 


1.3-2.3 


1.3-1.6 


ECAL miscalibration 


< 1 


< 1 


< 1 


< 1 


Total on F to t a i,r, 


4.4-13 


5.4-11 


17-49 


4.4- 12 



by removing from the measurement all tracks found within the same event with similar 
momentum vectors. It is observed that the charged energy flow drops by less than 1% 
for all i] bins and event classes in case of both data and simulation. For the final results, 
a global 1% systematic uncertainty for multiply reconstructed tracks is conservatively 
assigned. The effect of spurious tracks is estimated in simulation by determining the 
energy flow carried by reconstructed tracks which cannot be associated with particles 
at generator level and conservatively assuming that the simulation describes the rate of 
spurious tracks with a precision of 20%. The resulting systematic uncertainties are found 
to vary between 0.2% and 2%. 

It has been checked that reversing the LHCb magnet polarity has only an influence on 
the measurements of the charged energy flow for the diffractive enriched event class, which 
mainly consists of low-multiplicity events. Here, the corresponding effect is assigned as a 
systematic uncertainty. 

Events with more than one reconstructed primary vertex are vetoed in the analysis in 
order to suppress pile- up contamination. Its residual effect is estimated to be 1.7% by tak- 
ing the efficiencies to accept pile-up events from simulation. This factor is included in the 
normalisation of the energy flow and conservatively taken as the systematic uncertainty. 

The total energy flow acquires an additional uncertainty from the variation of the 
R g en,r) and k v factors between the Pythia tunes and the extrapolation procedure used for 
the k v factors in two highest r\ bins. No systematic uncertainty is assigned to account for 
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Table 2: Charged energy flow for all event classes and rj bins with the corresponding systematic 
uncertainties. The statistical uncertainties are insignificant and not listed. All values are in 
GeV per unit pseudorapidity interval. 





J-11L<1 LLol V \Z, 


r4 a T*H 
±±ObL H 


T 11 TTT C 1 ^ 1 ATP 


IxToTi _H 1TTT n r*"i" i"\rf^ 


range 


minbias 


scattering 


enriched 


enriched 


1.9 < 7] < 2.2 


12±1 


38±4 


4±2 


13±1 


2.2 < T] < 2.5 


16±1 


51±4 


5±2 


18±1 


2.5 < 7] < 2.8 


21±1 


66±4 


7±2 


23±1 


2.8 < 7] < 3.1 


28±1 


84±5 


9±3 


29±1 


3.1 < r] < 3.4 


36±2 


107±6 


12±3 


38±2 


3.4 < r] < 3.7 


46±2 


134±7 


17±4 


49±2 


3.7 < r] < 4.0 


59±2 


163±8 


23±5 


62±2 


4.0 < r] < 4.3 


74±3 


198±10 


32±7 


79±3 


4.3 < T] < 4.6 


90±4 


223±12 


42±7 


95±4 


4.6 < rj < 4.9 


114±6 


260±13 


58±10 


120±6 



Table 3: Total energy flow for all event classes and r\ bins with the corresponding systematic 
uncertainties. The statistical uncertainties are insignificant and not listed. All values are in 
GeV per unit pseudorapidity interval. 



Pseudorapidity 


Inclusive 


Hard 


Diffractive 


Non-diffr active 


range 


minbias 


scattering 


enriched 


enriched 


1.9 < 7] < 2.2 


18±2 


56±6 


4±2 


20±2 


2.2 < 7] < 2.5 


26±2 


79±7 


6±2 


28±2 


2.5 < 7] < 2.8 


36±2 


104±7 


10±3 


39±2 


2.8 < 7] < 3.1 


49±3 


135±9 


16±5 


52±3 


3.1 < 7] < 3.4 


61±3 


167±9 


20±5 


65±3 


3.4 < 7] < 3.7 


76±3 


207±11 


28±7 


81±4 


3.7 < 7] < 4.0 


96±4 


250±15 


38±9 


102±4 


4.0 < 7] < 4.3 


120±5 


301±18 


50±11 


127±6 


4.3 < rj < 4.6 


146±7 


334±20 


66±12 


154±7 


4.6 < 7] < 4.9 


186±10 


386±22 


91±15 


194±10 



inaccuracies of the Lund string model in describing the ratio of the neutral energy flow to 
the charged energy flow. The uncertainties associated with the ECAL energy calibration 
and photon reconstruction efficiency also affect the accuracy of the total energy flow. To 
account for a possible difference in the photon reconstruction efficiency between simulation 
and data, a global 3.7% systematic uncertainty is assigned following the analysis presented 



in Ref. 35 . The energy calibration of the ECAL has been studied by measuring the 
invariant masses of diphoton resonances (n° — > 77 and 77 —> 77) and has been assigned 
a global systematic uncertainty of 1.5%. Both uncertainties are scaled with a factor 
Fneuttf/Ftot^r, and are listed in Table [l] 
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Other potential sources of systematic uncertainties such as momentum- and rj- 
smearing, effect of the beam crossing angle, neglecting the masses of charged particles, 
pollution from elastic scattering and beam-gas interactions have been studied as well. 
Their impacts on the accuracy of the measurements are found to be negligible. 

The total systematic uncertainties on the corrected charged and total energy flow are 
listed in Tables [2] and [3] for all event classes and rj bins. It should be noted that the 
uncertainties are strongly correlated between the bins. 

5 Results 

The fully-corrected measurements for the charged energy flow are shown in Fig. [I] for each 
event class together with the generator level predictions given by the Pythia tunes and 
the corresponding systematic uncertainties. By comparing experimental results obtained 
for different event classes one can clearly see that the amount of energy flow strongly 
correlates with the momentum transfer in an underlying pp inelastic interaction. The 
charged energy flow rises more steeply with pseudorapidity in data than predicted by the 
majority of the Pythia tunes. As a consequence, the discrepancy between the measure- 
ments and generator level predictions increases towards large 77 rising to 20% in the last 
r\ bin. At lower rj the data are reasonably well described by the Pythia tunes. This is 
the case for all event classes except the diffractive enriched one. For the latter, there is a 
rather large inconsistency between the data and MC models for lower values of rj, while at 
large r] the measurements are well described by the Pythia 8.135 generator with default 
parameters. However, this model overestimates the charged energy flow in the case of 
hard scattering events over the entire pseudorapidity range of the measurements. 

Figure [2] illustrates the charged energy flow along with the predictions given by the 
cosmic-ray interaction models. It is interesting to note that the measurements performed 
with inclusive minimum-bias and non-diffractive enriched events are well described by the 
Epos 1.99 and Sibyll 2.1 models, while the QgsjetOI and QgsjetII-03 generators 
overestimate the charged energy flow for these event classes. The latter also occurs at 
large r] in the case of hard scattering interactions for all cosmic-ray interaction models 
except the QgsjetII-03. The diffractive enriched charged energy flow is underestimated 
by all cosmic-ray interaction models. 

The total energy flow is shown for each event class in Fig. [3] along with the generator 
level predictions given by the Pythia tunes and the corresponding systematic uncer- 
tainties. It can be clearly seen that all Pythia 6 tunes underestimate the amount of 
energy flow at large pseudorapidity for all event classes. The Pythia 8.135 generator 
gives the best description of the measurements performed with inclusive minimum-bias, 
diffractive and non-diffractive enriched events among the Pythia tunes, except for the 
pseudorapidity range 1.9 < 77 < 2.5. None of these models provide an accurate description 
of the energy flow measured with hard scattering events. The predictions given by the 
LHCb and Perugia NOCR tunes for non-diffractive enriched and hard scattering events 
are rather similar, while the Perugia tune significantly underestimates the energy flow 
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Figure 1: Charged energy flow as a function of i] for all event classes as indicated in the figures. 
The corrected measurements are given by points with error bars, while the predictions by the 
Pythia tunes are shown as histograms. The error bars represent the systematic uncertainties, 
which are highly correlated between the bins. The statistical uncertainties are negligible. The 
ratios of MC predictions to data are shown in addition. 



10 



> 
o 



LU 



<D 

CD 
i_ 
cG 

O 



05 
CO 

Q 



180 
160 
140 
,120 
100 
80 
60 
40 
20 




— • — Data 

EPOS 1.99 

QGSJET01 

QGSJETII-03 

SYBILL 2.1 
inclusive minbias events 



^•350 
a> 

±±300 



LU 



250- 



LHCb 

\fs=7 TeV 



f\ 200 
&150 

i 

CO 

O 100 



50 




— • — Data 

EPOS 1.99 

QGSJET01 

QGSJETII-03 

SYBILL 2.1 
hard scattering events 



LHCb 
{s=7 TeV 

I i i i i I 



1.4 




I Systematic Uncertainty _ 


1.2 
1 






0.8 






0.6 


i i i i i i 


I i i i i I i i i i I i i i i I i i i 







m 


1.4 


a 




o 


1.2 








1 




0.8 




0.6 



Systematic Uncertainty 



2.5 



3.5 



4.5 



2.5 



3.5 



4.5 



> 70 FY 

$ - — •— Data 

" 60- EPOS 1.99 

-g" : QGSJET01 

Uj QGSJETII-03 

SYBILL 2.1 
diffractive enriched events 



co 30 — 




i i i i | i i i i | i i i i | i 
— • — Data 

EPOS 1 .99 

QGSJET01 

QGSJETII-03 

SYBILL 2.1 
non-diffractive enriched events 




LHCb 
^=7 TeV 



_L 



0.5- 



1.4 




Systematic Uncertainty _ 


1.2 




1 






0.8 






0.6 


i i i i I i i i i 


i i i i i i I i i i i i i 



2.5 



3.5 



4.5 



2.5 



3.5 



4.5 



Figure 2: Charged energy flow as a function of r/ for all event classes as indicated in the figures. 
The corrected measurements are given by points with error bars, while the predictions by the 
cosmic-ray interaction models are shown as histograms. The error bars represent the systematic 
uncertainties, which are highly correlated between the bins. The statistical uncertainties are 
negligible. The ratios of MC predictions to data are shown in addition. 



11 



^250 

CD 

o 

€200- 

UJ 

■o 

Z S 1 50 



100 



50- 



— • — Data 

PYTHIA6 LHCb 

PYTHIA6 Perugia 

PYTHIA6 Perugia NOCR 

- PYTHIA 8.1 35 default 
inclusive minbias events 




LHCb 

\fs=7 TeV 



05 




CO 


1.4 


Q 




O 


1.2 








1 




0.8 




0.6 



I Systematic Uncertainty 



2.5 



3.5 



4.5 



> 

$100- 

€ : 

W 80- 

"D 



— • — Data 

PYTHIA6 LHCb 

- PYTHIA 8.130 diffraotive 
PYTHIA 8.135 default 
diffractive enriched events 



^ 60- 
$40- 



20- 





Q 
O 



US 



4= 



I 



I 



I" 



r — ' LHCb 
is=7 TeV 

I i i i i I 



Systematic Uncertainty 



0.5- 



2.5 



3.5 



4.5 



^•450 

CD 

O 400 



$ 350 - 
UJ 

~°-300 

z 

^ 250 

1200 
F 

150 
100 
50 




— • — Data 

PYTHIA6 LHCb 

PYTHIA6 Perugia 

PYTHIA6 Perugia NOCR 

- PYTHIA 8.135 default 
hard scattering events 



LHCb 
{s=7 TeV 

I i i i i I 







<3 


r.4 


Q 




O 


f.2 








/ 




0.8 




0.6 



Systematic Uncertainty 



2.5 



3.5 



4.5 



^■250r 
CD 

O 



T 



— • — Data 

PYTHIA6 LHCb 

^200[ PYTHIA6 Perugia 

[Q - PYTHIA6 Perugia NOCR 

"D I - - PYTHIA 8.135 default 

non-diffractive enriched events 



~150- 



700- 



50- 



.±r r 



LHCb 
^=7 TeV 



B 




(0 


1.4 


Q 




O 


1.2 








1 




0.8 




0.6 



11 111 






Systematic Uncertainty _ 








i=N i 


i 


I I 


1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 



2.5 



3.5 



4.5 



Figure 3: Total energy flow as a function of rj for all event classes as indicated in the figures. 
The corrected measurements are given by points with error bars, while the predictions by the 
Pythia tunes are shown as histograms. The data obtained with extrapolated factors are 
shown in grey. The error bars represent the systematic uncertainties, which are highly correlated 
between the bins. The statistical uncertainties are negligible. The ratios of MC predictions to 
data are shown in addition. 
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Figure 4: Total energy flow as a function of rj for all event classes as indicated in the figures. 
The corrected measurements are given by points with error bars, while the predictions by the 
cosmic-ray interaction models are shown as histograms. The data obtained with extrapolated 
k v factors are shown in grey. The error bars represent the systematic uncertainties, which are 
highly correlated between the bins. The statistical uncertainties are negligible. The ratios of 
MC predictions to data are shown in addition. 
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for all event classes. For diffractive enriched events, the inconsistency between the data 
and the prediction given by the LHCb tune is found to be rather large throughout the 
entire pseudorapidity range 1.9 < 77 < 4.9, while the Pythia 8.135 generator with default 
parameters gives a good description of the corresponding energy flow at large pseudora- 
pidities. 

Figure [4] illustrates the total energy flow together with the predictions given by the 
cosmic-ray interaction models. It is observed that the Sibyll 2.1 generator gives the best 
description of the energy flow measured with inclusive minimum-bias and non- diffractive 
enriched events at large 77. The predictions given by the Epos 1.99 generator for these 
event classes also describe the measurements reasonably well. In the case of hard scatter- 
ing interactions, the best description of the data at large rj is given by the QgsjetII-03 
generator. The total energy flow measured with diffractive enriched events is underesti- 
mated at large 77 by all cosmic-ray interaction models used in this study. The measure- 
ments of the charged and total energy flow are summarised in Tables [2] and [3] for all event 
classes and 77 bins. 

The results obtained in this study cannot be directly compared with the measurements 
performed by the CMS collaboration [5] , since different event selection criteria are applied 
in the analyses. Nevertheless, both measurements demonstrate that the energy flow is 
underestimated by Pythia 6 tunes at large pseudorapidity. 

6 Conclusions 

The energy flow is measured in the pseudorapidity range 1.9 < 77 < 4.9 with data col- 
lected by the LHCb experiment in pp collisions at a/s = 7 TeV for inclusive minimum-bias 
interactions, hard scattering processes and events with enhanced or suppressed diffractive 
contribution. The primary measurement is the energy flow carried by charged particles. 
For the measurement of the total energy flow, a data-constrained MC estimate of the neu- 
tral component is used. The energy flow is found to increase with the momentum transfer 
in an underlying pp inelastic interaction. The evolution of the energy flow as a function 
of pseudorapidity is reasonably well reproduced by the MC models. Nevertheless, the 
majority of the Pythia tunes underestimate the measurements at large pseudorapidity, 
while most of the cosmic-ray interaction models overestimate them, except for diffractive 
enriched interactions. For inclusive and non-diffractive enriched events, the best descrip- 
tion of the data at large 77 is given by the Sibyll 2.1 and Pythia 8.135 generators. 
The latter also provides a good description of the energy flow measured with diffractive 
enriched events, except for the pseudorapidity range 1.9 < 77 < 2.5. The comparison 
shows that the absence of hard diffractive processes moderates the amount of the for- 
ward energy flow meaning that their inclusion is vital for a more precise description of 
partonic interactions. It also demonstrates that higher-order QCD effects as contained in 
the Pomeron phenomenology play an important role in the forward region. None of the 
event generators used in this analysis are able to describe the energy flow measurements 
for all event classes that have been studied. 
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